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ABSTRACT 


Computer  graphics  is  particularly  appropriate  foi 
studying  the  time  response  of  two-dimensional  mechanical 
dynamic  systems.  This  paper  discusses  the  interactive 
use  of  a  DEC  338  graphics  terminal  for  input  and  display 
related  to  a  generalized  FORTRAN  IV  applications  program 
named  DAMN- - Dynamic  Analysis  of  Mechanical  Networks. 

The  DEC  338  is  remotely  timeshared  to  the  University  of 
Michigan  IBM  System  360/67  central  computer,  operating 
under  MTS  (Michigan  Terminal  System). 


1.  INTRODUCTION 


This  paper  deals  with  the  use  of  computer  graphics 
in  an  attempt  to  apply  computer-aided  design  to  realistic 
machinery.  Here  the  term  machinery  denotes  mechanical 
machinery,  including  such  devices  as  high-speed  card 
punches,  automotive  suspensions,  washing  machines, 
tractor  backhoes,  and  control  linkages.  Approximately 
one  tenth  of  all  engineers  are  concerned  primarily  with 
design  of  this  class  of  machinery.1  Characteristics  of 
realistic  machinery  include  the  following: 

Large-scale  multidimensional  movement.  In  con¬ 
trast  to  the  vibration  of  structural  systems, 
the  time  response  of  machinery  may  involve 
translational  movements  comparable  to  the  size 
of  the  machine  itself.  Angular  coordinates 
may  move  through  complete  rotations.  Typically 
these  displacements  occur  in  multiplanar  two- 
dimensional  space,  although  often  machinery  is 
designed  for  motion  which  can  only  be  described 

t h ree-dimensionally. 

\ 

2  Detailed  constraint.  Most  realistic  machinery 
does  not  conform  to  the  lumped  mass- spring- 
damper  model  used  to  explain  vibration.  Typi¬ 
cally  parts  are  interconnected,  not  only  by 
elastic  and  damping  elements,  but  by  other  parts 
as  well.  These  part-to-part  contacts  frequently 
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form  complicated  loop  structures.  Moreover 
the  contacts  themselves  are  varied.  Some  are 
essentially  simple  hinges  or  keyways  and  are 
fixed  in  location  relative  to  the  contacting 
parts.  However,  others  shift  location  through¬ 
out  the  motion  of  the  machine.  An  example  is 
the  contact  between  a  cam  and  a  follower. 

3.  Multifreedom  dynamics.  Despite  the  extensive 
constraint  present  in  machinery,  it  is  still 
necessary  to  regard  it  as  having  more  than 

zero  degrees  of  freedom.  (The  classical  approach 
to  machinery  involves  just  this  assumption- - that 
machinery  is  kinematic  in  nature,  with  zero  free¬ 
dom  once  all  input  motions  are  specified.)  Multi 
freedom  behavior  may  be  introduced  by  elastic 
elements  such  as  timing  belts  or  flexible  beams, 
by  motors  which  lack  the  power  to  bring  the 
machine  almost  instantly  up  to  speed,  or  by 
unusual  inteiaction  of  a  normally  rigid  system 
(automobile  crash) 

4.  Nonlinear  elements.  Typically  the  elements  of 
realistic  machinery  are  not  ideal  springs  and 
dampers  but  instead  such  items  as  bushings  or 
pneumatic  elements  For  large-scale  motion, 
these  cannot  be  represented  in  a  single  linear 


model 
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5  Intermi t tency .  Often  a  machine  will  change 
its  structural  identity,  perhaps  irregularly, 
at  times  during  its  motion.  Higher  pair  contacts 
may  open  because  of  inertial  loads,  or  parts 
which  are  normally  fixed  in  position  may  begin 
to  move  due  to  overload. 

6.  Multidisciplinary  designs.  Mechanical  designers 
tend  to  actively  couple  electronic,  optical, 
fluid,  and  other  diverse  elements  in  their 
machine  design. 

7.  Additional.  Other  deviant  characteristics  of 
machinery  will  typically  be  observed  according 

to  the  product  with  which  a  machine  is  associated. 

For  example,  thread  manufacturing  machinery  has 
parts  which,  in  effect,  change  mass  with  time 
as  more  thread  is  wound  onto  a  spool. 

Effective  computer-aided  design  of  any  type  of  engin¬ 
eering  system  requires  at  its  foundation  a  generalized  pro¬ 
gram  which  will  accept  not  only  changes  in  the  values  of  the 
parameters  of  a  specific  device  but  will  also  accept  changes 
in  the  identity  of  the  device  itself.  As  suggested  by  the 
foregoing  characterization,  the  development  of  a  true 
generalized  program  for  computer-aided  design  of  realistic 
machinery  is  an  imposing  task.  In  fact,  the  program  dis¬ 
cussed  in  this  paper  is  somewhat  more  limited  than  the  ideal. 

It  accommodates  large-scale  motion  but  only  for  two-dimensional 
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systems.  Ihis,  however,  sets  it  apart  from  most  dynamics 
programs  designed  for  the  study  of  vibration.  It  accommo¬ 
dates  the  constraint  of  part  loop  closure  and  therefore 

applies  to  two-dimensional  kinematic  systems  such  as  those 
2 

covered  by  KAM  and  the  generalized  matrix  kinematic  pro- 

3 

grams  written  at  Northwestern  University.  KAM  and  the 
Northwestern  program  cover  three-dimensional  systems  also, 
but  do  not  extend  to  dynamics.  Because  of  problems  of 
detail,  the  present  program  is  limited  in  covering  higher 
pairs  and  nonlinear  elements,  and  has  reasonable  limits  on 
problem  size.  But  it  does  provide  a  working  basis  for  a 
large  fraction  of  computer-aided  machine  design.  The  pro¬ 
gram  algorithms  were  aimed  toward  a  system  representation 
for  machinery  similar  to  that  which  has  been  so  successful 
for  electrical  and  structural  systems.  Because  of  the 
complexity  previously  cited,  a  complete,  consistent  identi- 

4 

fication  with  either  network  state  variable  techniques  or 
bond  graph^  techniques  has  not  yet  been  made.  However  the 
identification  is  sufficiently  close  that  the  program  has 
been  named  DAMN  (Dynamic  Analysis  of  Mechanical  Networks). 

Computer  graphics  is  a  most  helpful  adjunct  in  applying 
the  DAMN  program  to  machine  design.  Graphic  output  allows 
the  designer  a  very  rapid  appreciation  of  the  relation  of 
the  several  problem  coordinates  as  the>  perform  in  two- 
dimensional  space  He  can  immediately  sense  the  phase  re¬ 
lations,  accelerations,  possible  interference  problems, 
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and  other  features  which  are  obviously  difficult  to  inter¬ 
pret  from  a  collection  of  individual  plots  of  the  variables. 
Because  of  the  effectiveness  of  graphic  output  there  is  an 
increased  emphasis  on  the  worth  of  an  interactive  capability 
to  the  designer.  The  program  has  been  designed  to  allow 
immediate  change  in  either  the  parameters  defining  the  de¬ 
sign  of  a  specific  system  (masses,  stiffness  coefficients, 
geometric  lengths,  etc.)  or  the  topological  data  which  dis¬ 
tinguish  one  system  from  another. 

Figure  1  shows  generally  the  relation  between  program, 
remote  graphic  input/output ,  the  University  of  Michigan 
central  computer  facility,  and  the  DAMN  program  itself. 

The  DEC  338  display  terminal  of  course  includes  a  PDP-8 
computer  connected  to  a  graphic  display  unit.  This  is 
connected  in  full-duplex  mode  over  normal  telephone  lines 
to  another  computer,  also  a  PDP-8,  modified  for  communica¬ 
tions  use  and  termed  the  Data  Concentrator.  In  effect, 
the  graphic  terminal  can  communicate  with  the  DAMN  program 
at  any  level  for  flexible  simulation  of  a  specific  device 
or  set  of  devices.  The  DAMN  program  itself  is  written  en¬ 
tirely  in  FORTRAN  IV  It  is  resident  on  disc  storage  at 
the  University  of  Michigan  Computing  Center  and  is  executed 
on  the  University's  IBM  360/67  computing  system**.  The 
graphics  output  program  also  shown  is  written  primarily  in 
FORTRAN  IV,  but  is  importantly  dependent  on  the  DF  routines 

if 

developed  under  the  Concomp  Project  for  the  support  of  the 
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remote  graphic  terminal.10  At  present,  light-pen  input  is 
not  used  in  a  major  way.  Instead,  desired  parameter  changes 
are  typed  in  from  a  teletypewriter  unit  on-line  to  the  DEC 
338  display  terminal. 

The  DAMN  program  was  created  as  part  of  a  project 
sponsored  primarily  by  the  Ford  Motor  Company  for  the 
computer  analysis  of  mechanical  dynamic  systems.  The 
project  was  substantially  assisted  also  by  the  Concomp 
Project,  itself  sponsored  by  the  Advanced  Research  Projects 
Agency  of  the  U.S.  Department  of  Defense.  The  DAMN  program 
is  viewed  as  a  prototype  program  useful  in  giving  students 
an  appreciation  of  computer-aided  design  of  mechanical 
systems  and  in  revealing  to  them  the  worth  of  analytical 
methods  and  computer  characteristics  for  subsequent  pro¬ 
gramming  of  a  more  advanced,  comprehensive  system.  How¬ 
ever,  experience  with  DAMN  to  date  has  been  very  favorable 
particularly  the  level  of  detail  and  speed  with  which  the 
response  of  mechanical  systems  can  be  determined,  and  the 
characteristics  of  graphic  display  of  mechanical  systems. 

In  Section  2  the  analytical  background  and  program 
structure  of  DAMN  will  be  reviewed  more  specifically  Sec¬ 
tion  3  will  detail  the  graphic  programming  and  system 
facilities  and  discuss  experience  with  specific  mechanical 
systems  which  have  been  simulated  in  graphic  computer  runs. 
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2.  BACKGROUND  OF  PROGRAM 

2 . 1  Cons t ituitive  Elements 

The  DAMN  program  is  designed  to  represent  machinery 
systems  by  a  simultaneous  set  of  second-order  ordinary 
differential  equations.  The  essential  programming  problem 
is  the  development  of  the  coefficient  terms  in  these  equa¬ 
tions  regardless  of  the  machinery  system  that  has  been  de¬ 
fined  to  the  program  via  the  input.  To  facilitate  this,  a 
lumped  parameter  view  of  machinery  is  taken  in  which 
machinery  is  seen  to  consist  of  eight  constituitive  ele¬ 
ments  : 

(1)  Parts .  Rigid,  inextensible  objects  of  zero  or 
finite  size,  mass,  and  moment  of  inertia* 

(2)  Fields  Conservative  or  nonconservative  forces 
applied  between  one  part  and  another. 

(31  Generators .  Motion  sources,  inputting  either  a 
rotational  or  translational  motion. 

(4)  Rotations  Contacts  allowing  a  single  degree 
of  freedom  of  rotation  of  one  part  relative  to 
another 

(5)  Translations  Contacts  allowing  a  single  degree 
of  freedom  of  translation  of  one  part  relative 
to  another. 

(6)  Markers  Points,  unit  vectors,  or  the  combination 
of  a  point  and  unit  vector,  fixed  trans lational ly 
or  in  orientation  relative  to  a  specified  part. 
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(7)  Keyed  Assemblies.  Assemblies  of  one  or  more 
parts  related  entirely  by  translational  contact 

(8)  Loops  Seq"ences  of  parts  and  contacts  forming 
an  independent  closed  chain 

All  data  input  to  the  DAMN  program  is  categorized  under 
the  first  six  elements.  Categorization  under  the  elements 
"Keyed  Assembly"  and  "Loop"  is  made  internally  by  the  pro¬ 
gram,  but  the  user  is  not  expected  to  recognize  these  ele¬ 
ments  themselves.  The  user  also  specifies  the  output  he 
desires  in  terms  of  the  relative  displacement,  velocity, 
and  acceleration  of  selected  pairs  of  markers.  The  program 
is  being  developed  also  to  compute  reaction  force  at  rota¬ 
tions  and  translations 

Based  on  the  data  input  in  the  preceding  categories, 
the  program  constructs  a  set  of  second-order  ordinary 
differential  equations  and  numerically  integrates  these 
through  the  range  of  time  specified  by  the  user  In 
general,  numerical  integration  is  a  slower  process  than 
the  direct,  algebraic  computation  ordinarily  possible  for 
kinematic  systems,  but  it  is  a  practical  necessity  for 
large  nonlinear  systems  Most  of  the  equations  con¬ 
structed  correspond  to  Lagrange  equations  One  of  these 
is  constructed  for  the  angle  of  each  keyed  assembly  in 
the  system,  and  one  for  the  displacement  of  each  trans¬ 
lation  Thus,  the  coordinates  for  two-dimensional 
mechanical  systems  are  chosen  as  keyed  assembly  angles 
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and  translational  displacements.  Additional  differential 
equations  are  created  to  represent  constraint. 


2 . 2  Definition  of  Constraint 

Figure  2.1  shows  a  revealing  though  very  academic 
example  of  a  machine.  The  coordinates  for  this  system 
are  the  angles  of  the  two  movable  keyed  assemblies  2  and 
3,  plus  the  three  translational  coordinates  Sjj,  s2i* 
and  s^.  For  any  situation  of  forces  acting  between 
the  parts  in  this  example,  the  system  is  dynamic  with 
three  degrees  of  freedom.  However,  the  fact  that  the 
parts  are  connected  in  a  closed  chain  or  loop  means 
that  the  coordinates  are  interrelated  by  the  following 
equation . 


r  *  s  ♦  r  ♦  r 
r2  1  b  2 1  r2  2  31 


-+ 
♦  s 


3rr32  +  rll  +  Sll+r12 


(1) 


Equation  (1)  simply  states  that  the  sum  of  part  vectors  r 
and  translation  vectors  s  around  the  closed  loop  is  0 

A  /V 

The  components  of  this  equation  in  the  directions  i  and  j 
correspond  to  two  scalar  equations  in  the  coordinates  ©2, 
0,,  S2j,  S-,  ,  and  s^.  When  the  problem  of  Figure  2.1 
is  represented  by  DAMN,  the  second  derivative  of  the  com¬ 
ponents  of  Eq  (l)  is  included  as  a  pair  of  constraint 
differential  equations  supplementing  the  Lagrange  equations 


As  a  limiting  case  the  constraint  conditions  are  the 
conditions  normally  used  to  solve  kinematic  problems  (zero 


translations  are  denoted  more  conveniently  by  small  triangles 
The  angles  02  and  0  3  shown  are  the  angles  of  the  keyed  assem¬ 
blies  relative  to  ground. 
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degree  of  freedom).  For  example,  in  Figure  2.1,  if  the 
slips  s21  and  s31  are  prevented  and  the  angle  of  keyed 
assembly  number  2  is  specified  as  a  function  of  time, 
the  system  becomes  a  kinematic  slider  crank  mechanism, 
as  shown  in  Figure  2.2.  In  this  case  only  the  con¬ 
straint  equations  are  necessary  to  solve  the  system, 
either  by  direct  evaluation  of  algebraic  equations 
as  F.q.  (1),  or  by  numerical  integration  of  a  derivative 
of  Eq  (1)  . 

For  a  program  to  simulate  three-dimensional  machinery 
another  basic  constraint  condition  must  be  included-- 
specifically  the  condition  of  angular  constraint  around 
a  closed  loop  of  keyed  assemblies. 

2 . 3  Dynamic  Condi t ions -- Lagrange 1 s  Equation 

The  dynamic  condition  on  which  the  DAMN  program  is 
based  is  the  Principle  of  Least  Action,  or  Hamilton's 
Principle  As  applied  to  mechanical  systems  this  states 
that  a  system  will  respond  so  as  to  minimize  the  integral 
over  time,  of  the  work  due  to  all  effects,  by  the  surround¬ 
ings  on  the  system.  In  equation  form  this  condition  is 
presented  as  the  variation  of  an  integral  equal  to  0: 

6  /0tw(^i*qi't)  =  0  (2) 

By  use  of  variational  mathematics  this  condition  is  de¬ 
veloped  to  the  following  set  of  differential  equations,  one 


seit  (.tq.  C1JJ  is  sufficient  for  determination  of  system 
time  response,  given  an  input  of  02  versus  time. 
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Here,tt.he  work  W  done  by  the  surroundings  on  the  system 
includes  kinetic  and  potential  energy  effects,  plus 
effects  due  to  nonconservative  forces  such  as  those 


of  damping  Ordinarily  these  are  the  only  effects 
considered  However,  a  basic  assumption  of  Eq.  (3) 
is  that  the  coordinates  q^  are  independent.  Because 
of  the  loop  constraints  mentioned  in  Section  2.2  this 
condition  is  not  normally  met  by  the  coordinates  recog¬ 
nized  by  DAMN  for  machinery  systems.  To  overcome  this 
difficulty  the  loops  are  treated  as  broken  for  purposes 
of  Lagrange's  equation.  Formally,  the  coordinates  qi 
are  therefore  independent  as  required  However,  the 
Lagrange  equations  thereby  produced  are  accompanied  by 
second-order  loop  constraint  equations  which  as  a  practical 
matter  insure  that  the  loops  remain  closed.  Also  the 
work  terms  in  the  Lagrange  equations  now  include  Lagrange 
multiplier  terms  for  each  component  of  loop  constraint. 


These  have  the  form  A .  4> .  where  the  X.  are  reaction  forces 

J  J  J 


and  the  <t>  .  are  loop  constraint  functions. 


The  ultimate 


form  of  Lagrange's  equation  with  constraint  is  the  following: 
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As  it  stands  Eq.  (4)  is  difficult  to  program  for 
general  application.  What  is  instead  required  is  a 
set  of  relations  directly  of  the  form  of  second-order 
ordinary  differential  equations.  The  principal  problem 
here  is  the  complicated  set  of  operations  performed  on 
the  kinetic  energy  terms.  By  viewing  a  mechanical  sys¬ 
tem  in  terms  of  its  constituitive  elements  it  is  possible 
to  determine  general  expressions  for  the  Lagrange  equation 
kinetic  energy  terms,  for  the  mass  of  one  specific  part, 
part  m.  These  expressions  are  the  following: 
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These  relations  are  especially  interesting  when  the 
resemblance  is  noted  to  the  general  expression  for 
the  second  derivative  of  a  single  vector,  r. 

-►  •  •  A  -f  •  ••  /V  |  |  A  A 

r  =  0(k*r)  -  6  r  +  r  r  +  2r0(k*r)  (7) 

In  Eqs.  (5)  and  (6)  the  summation  over  i  is  a  summation 
over  all  the  keyed  assemblies  in  a  chain  from  ground  to 
the  part  whose  mass  is  be^.ng  considered.  The  summations 
over  j  are  summations  of  translational  terms  within  a 
given  keyed  assembly  i,  included  in  the  chain  mentioned. 
Figure  2.3  illustrates  the  elements  involved  in  this 
summation.  The  vectors  are  position  vectors  which 
span  the  ith  keyed  assembly,  from  the  contact  at  which 
the  chain  enters  to  the  contact  at  which  it  leaves. 

Two  equations  are  required,  (5)  and  (6),  because  the 
Lagrange  equation  must  be  written  both  for  coordinates 
which  are  angles  of  keyed  assemblies  0^  and  coordinates 
which  are  displacements  of  translational  contacts,  s ^ . 

As  far  as  the  assembly  of  inertial  terms  is  con¬ 
cerned,  at  each  integrating  step  DAMN  polls  each  part 
having  finite  mass  and  accumulates  coefficient  terms 
as  defined  by  Eqs.  (5)  and  (6). 

The  other  terms  in  Eq.  (4)  are  much  easier  to  evalu¬ 
ate.  The  Lagrange  multiplier  terms  Aj  4> ^  result  in  a 
matrix  of  coefficients  which  is  exactly  the  transpose  of 


the  matrix  of  coefficients  obtained  in  a  second-order 


Figure  2.3.  A  chain  of  keyed  assemblies  from  ground  to  the 
center  of  mass  of  a  particular  part.  In  constructing  the 
Lagrange  equation  kinetic  energy  terms  for  each  part  m  (Eqs. 
(5)  and  (6)),  summations  are  made  over  the  chain  of 
keyed  assemblies  to  the  part  center  of  maSs,  and  over  the  in¬ 
dividual  translations  within  each  keyed  assembly  in  the  chain. 


constraint  condition.  Forces  due  to  potential  energy 
effects  (operations  on  V)  are  simply  included  in  the 
collection  of  nonconservative  force  terms  Q^.  Table  1 
shows  the  construction  of  the  force  terms,  differing 
according  to  whether  the  Lagrange  equation  coordinate 
is  angular  or  translational  and  whether  the  actual 
applied  force  is  a  translational  force  or  a  rotational 
force  (torque) . 

The  forces  included  in  the  term  are  distinquished 
as  being  displacement-  or  velocity-dependent,  not  acceler¬ 
ation-dependent.  This  term  therefore  includes  spring 
forces,  damping  forces,  and  any  other  force  whether  linearly 
or  nonlinearly  dependent  on  displacement  and  velocity.  One 
standard  force-type  available  in  DAMN  is  a  "wall  force"  which 
has  magnitude  zero  until  a  specified  pair  of  points,  thought 
of  as  the  terminals  of  a  force  field,  move  within  a  dis¬ 
tance  comparable  to  an  integration  step  of  each  other.  The 
force  function  then  increases  very  rapidly  for  very  small 
decrease  in  point  separation.  In  effect  the  points  are 
prevented  from  exact  coincidence.  This  particular  force 
type  provides  a  basis  for  including  a  wide  variety  of 
intermittent  mechanisms  and  impact  situations  within  DAMN. 

At  present  DAMN  is  not  organized  to  accept  a  complete 
variety  of  one-of-a-kind  forces.  Such  acceptance  requires 
allowing  the  user  to  input  force  functions  in  either  sub¬ 
routine  or  tabulated  form,  certainly  a  feasible  extension 


TABLK  1 


Form  of  the  Generalized  Force  Terms 
in  the  Lagrange  Equation 


Type  of 

App 1 i ed 

Force 

Di sp 1 acement 

Coord i nate  ,  s  . 

Up 

Angular 

Coordinate,  0^ 

Trans  1  at lona 1 

Force 

-* 

Fk 

-*  A 

Fk  s  z.p 

Fk  (k  X  b£) 

Rotat l ona 1 

Force  (torque) 

-> 

Tk 

Zero 

BB 

F^  =  kth  actual  applied  translational  force  in  the 
system 

1  j.  =  kth  actual  applied  rotational  force  (torque)  in 
the  system,  applied  between  a  pair  of  parts  1 
and  j 

=  angle  of  the  Jet h  keyed  assembly 
Snp  ‘  the  pth  translation  within  the  y,th  keyed  assembly 
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of  the  program.  A  similar  restriction  now  exists  for 
one-of-a-kind  motion  for  generators  and  part  contacts 
dependent  on  surface  profile.  When  the  reaction  force 
computation  capability  is  completed  DAMN  will  accept 
applied  forces  which  are  themselves  dependent  on 
reaction  force  (Coulomb  friction  force).  Meanwhile 
this  force  type  is  excluded. 


The  essential  form  of  the  dynamic  and  constraint 
conditions  necessary  to  DAMN  can  be  represented  as  an 
Cq.  (8): 
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Here  the  differential  equations  are  represented  as  a  set 
of  simultaneous  linear  algebraic  equations  in  the  variables 
q  ^  and  a  j  The  coefficient,  matrix  is  partitioned  in  four 
submatrices  The  upper  left  submatrix  contains  the  co¬ 
efficients  A  of  the  second-order  acceleration  terms 
from  Lagrange's  equation  Specifically  these  are  radial 
and  tangential  acceleration  coefficients  The  main 
diagonal  coefficients  of  this  matrix  are  constants,  and 
the  submatrix  is  symmetric  However,  the  off-diagonal 
terms  are  nonlinearly  dependent  on  problem  coordinates, 
specifically  the  angles  of  keyed  assemblies  The  upper 
right  and  lower  left  submatrices  happen  to  be  the  transpose 
of  each  other,  although  the  lower  left  submatnx  is  formed 
directly  from  the  second-order  terms  of  the  constraint 
conditions,  whereas  the  upper  right  submatrix  is  com¬ 
posed  of  Lagrange,  multiplier  coefficients  from  the 
lagrange  equations  themselves  These  coefficients 
are  also  nonlinearly  dependent  on  the  angles  of  the 
keyed  assemblies  The  column  vector  of  dependent  variables 
includes  the  second  derivatives  of  the  keyed  assembly 
angles  plus  the  second  derivatives  of  the  translational 
variables  The  lower  subvector  includes  the  Lagrange 
multiplier  terms  which  originate  from  Lagrange's  equation 
On  the  right  is  a  column  vector  of  elements  which  them¬ 
selves  are  large  summations  of  displacement-  and  velocity- 
dependent  forces,  applying  to  the  particular  Lagrange  or 
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constraint  equation  identified  with  their  row.  None  of 
these  terms  however  is  dependent  on  coordinate  second 
derivatives . 

The  actual  integration  process  takes  place  as  out¬ 
lined  in  the  lower  left  of  Figure  1.  First,  the  initial 
conditions  are  used  to  evaluate  all  of  the  coefficients 
in  Eq.  (8)  at  time  equal  zero.  Then  Eq.  (8)  is  solved 

as  a  simultaneous  set  of  linear  algebraic  equations  for 

••  •• 
the  q's  and  the  X's.  In  concept  the  q  terms  are  numeri- 

cally  integrated  to  determine  corresponding  velocity  and 
displacement,  then  fed  back  to  compute  the  coefficient 
terms  in  Eq  (8)  for  a  next  step.  At  present  the  inte¬ 
grating  program  used  is  HPCG  (Hamming  Predictor-Corrector 
Generalized)  available  from  the  IBM  Scientific  Subroutine 
Package  The  integration  takes  place  in  small  steps  with 
a  simultaneous  equation  solution  in  each  step  until  a 
use r- spec i fied  output  increment  is  completed.  At  that 
time  current  values  of  all  problem  coordinates  are  trans¬ 
mitted  to  the  graphics  output  program  and  processed  there 
for  immediate  display,  or  if  a  printed  output  has  been 
elected,  the  results  are  stored  on  disc  until  the  end 
of  the  entire  run,  at  which  time  they  are  printed  out. 

The  integration  process  outlined  here  is  simple, 
but  not  the  most  efficient  process  possible.  For  example, 
it  is  probably  feasible  to  integrate  without  having  to 
solve  Eq  (8)  at  every  integrating  step.  Also,  with  this 
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particular  set  of  variables  and  the  Lagrange  multiplier 
technique,  the  number  of  differential  equations  exceeds 
the  degree  of  freedom  of  the  system  by  the  amount  2c, 
where  c  is  the  number  of  scalar  constraints.  Alternate 
approaches  may,  in  fact,  yield  a  smaller  number  of  equa¬ 
tions  in  the  set  However,  the  Lagrange  equation  approach 
has  been  especially  easy  to  interpret  and  free  from  arti- 
fical  singularity  problems  because  of  the  close  physical 
correspondence  of  each  of  the  Lagrange  and  constraint 
equations  For  the  size  of  problems  considered  so  far 
(order  10  X  10),  integration  speeds  have  been  more  than 
adequate  for  analysis,  as  discussed  in  Section  3.  How- 
pver,  use  of  the  program  for  computer-controlled  syn¬ 
thesis  will  require  at  least  an  order  of  magnitude  in¬ 
crease  in  speed— probably  feasible 

2  4  Relation  of  Topological  and  Dimensional  Programs 

We  have  made  careful  distinction  between  topologi¬ 
cal  and  dimensional  subroutines  within  DAMN  Here,  the 
term  "topological  programs"  denotes  those  which  deter¬ 
mine  information  which  is  independent  of  part  dimensions, 
or  any  system  displacements  or  forces  A  good  example  of 
a  topological  program  is  the  subroutine  which  determines 
the  independent  loops  of  parts  and  contacts  in  an  arbitrary 
machine  system  The  results  of  this  program  are,  of  course, 
unaffected  by  the  state  of  the  machine  displacement  As 
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shown  in  Figure  1,  the  program  was  designed  so  that  all 
computations  which  are  topological  in  nature  are  per¬ 
formed  before  the  beginning  of  numerical  integration. 

The  results  of  these  computations  are  stored  as  integer 
information  in  a  set  of  program  tables,  termed  plex 
tables.  In  essence,  one  such  table  exists  for  each 
of  the  constituitive  elements  defined  in  Section  2.1. 

These  plex  tables  are  held  in  memory  throughout  both 
the  topological  and  dimensional  programs.  The  informa¬ 
tion  contained  considerably  increases  the  speed  with 
which  coefficient  values  can  be  determined  for  the  numeri¬ 
cal  integration  process. 

The  DAMN  program  is  network-br sed  in  the  sense  that 
most  of  the  topological  routines  use  network  mathematics 
corresponding  closely  to  that  used  in  electrical  or 
structural  analysis  problems  For  example,  the  minimum 
number  of  independent  part -contact  loops  is  determined 
by  processing  a  branch-node  incidence  matrix  where  parts 
are  interpreted  as  nodes  and  contacts  as  branches. 
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3.  GRAPHICS  SYSTEMS  AND  APPLICATIONS 

3 . 1  Systems  Viewpoint 

DAMN  is  a  set  of  FORTRAN  IV  subroutines  which  are 
executed  in  either  batch  or  timeshared  mode  on  the  Univer¬ 
sity's  IBM  360/67  central  computer.  All  such  programs  are 
under  the  control  of  a  supervisory  program  named  UMMPS 
(University  of  Michigan  Multiprogramming  System)  which 
supervises  a  timesharing  monitor,  MTS  (Michigan  Terminal 
System),  and  an  auxiliary  batch  monitor  HASP  (Houston 
Automatic  Spooling  Priority  System)  HASP  assigns 
priorities  to  batch  jobs.  In  general,  the  user  does 
not  have  access  to  UMMPS  or  HASP  and  must  communicate 

to  the  central  computer  through  MTS  by  means  of  the  MTS 

8 

command  language.  Figure  3.1  shows  the  interrelation  of 
these  programs  to  the  user.  Table  2  outlines  the  impor¬ 
tant  functions  of  UMMPS,  HASP,  and  MTS 

Currently  MTS  can  service  about  75  users  in  a  time- 
shared  mode  For  maximum  efficiency,  batch  jobs  are  first 
assigned  a  priority  by  HASP.  These  are  then  passed  to  MTS 
highest  priority  first,  for  execution  Batch  jobs  are  pro 
cessed  by  three  IBM  2540  card  read-punch  units,  one  IBM 
2501  card  reader,  and  five  IBM  1405  line  printers  Ter¬ 
minal  jobs  utilize  standard  IBM  and  AT&T  communications 
equipment  and  are  serviced  by  one  IBM  2702  and  one  IBM 
2703  communications  terminal. 


TERMINAL  BATCH 


ligure  3.1.  Relation  of  U  of  M  system  software  in  the  exe¬ 
cution  of  batch  and  timeshared  jobs  on  the  360/67  central 
computer.  Terminal  jobs  are  executed  according  to  MTS  com¬ 
mands.  Batch  jobs  are  first  assigned  a  priority  by  HASP, 
then  executed  in  timeshared  coordination  with  terminal  jobs 
Overall  control  is  accomplished  by  UMMPS. 


TABLE  2 


UMMPS 


MTS 


HASP 


Functions  of  the  Supervisory 


Programs  of 


The  University  of  Michigan 


IBM  360/67 


(University  of  Michigan  Multiprogramming  System). 
The  system  program  which  supervises  the  timesharing 
system.  Inaccessible  to  the  general  user.  Inter¬ 
faces  the  360  hardware  to  MTS.  Controls  the  load¬ 
ing  and  execution  of  user's  programs  and  I/O,  and 
processes  interrupts. 

(Michigan  Terminal  System).  User-oriented  system 
program.  Contains  assemblers,  translators,  and 
system  programs  which  the  user  accesses  by  issuing 
the  appropriate  command.  Monitors  both  terminal 
and  batch  jobs. 

(Houston  Automatic  Spooling  Priority  System). 
Interacts  with  UMMPS  to  assign  priority  to  batch 
jobs  for  effective  execution  of  programs.  Con¬ 
trol  is  eventually  passed  on  to  MTS  to  actually 
execute  program.  Inaccessible  to  general  user. 
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In  the  U  of  M  communications  network,  a  special 
control  unit  named  the  Data  Concentrator  can  be  accessed 
by  any  terminal  user.  The  Data  Concentrator  is  itself  a 
PDP-8  computer  containing  four  4096-word  banks  of  storage 
and  a  1 . 5 -mi c ro s econd  access  core  memory.  It  is  controlled 
by  a  supervisory  program  called  RAMP.  The  Data  Concentra¬ 
tor  was  developed  to  adapt  special  terminal  equipment  to 
the  timesharing  system,  and  its  applications  include  inter¬ 
facing  certain  non-IBM  special-purpose  communications 

9 

equipment  to  MTS.  For  DAMN  it  interfaces  a  DEC  338 
CRT  display  unit  to  the  central  computer,  making  possible 
the  graphic  display  feature.  Teletype  communication  to 
the  central  computer  through  the  Data  Concentrator  is 
essentially  the  same  as  through  the  IBM  equipment. 

3 . 2  The  Graphic  Terminal 

The  graphic  terminal  used  for  DAMN  output  is  a  DEC  338 
CRT  display  The  U  of  M  Concomp  Project  developed  the 
necessary  support  routines  to  interface  the  DEC  338  to 
the  central  computer.  A  user  communicates  with  the  central 
computer  over  voice-grade  telephone  lines  by  means  of  an 
ATGT  201  data  phone.  This  configuration  provides  the  user 
with  some  of  the  speed  advantages  of  the  central  computer, 
coupled  with  the  medium-cost  graphics  ability  of  the  338. 

The  338  also  serves  as  a  remote  graphics  terminal  time- 
shared  to  MTS  at  the  Ford  Motor  Company  in  Dearborn,  Michigan 
(25  miles  distant),  and  the  Whirlpool  Corporation  in  Benton 
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Harbor,  Michigan  (160  miles  distant).  These  terminals 
have  the  same  capabilities  as  those  located  at  the  Univer¬ 
sity  and  provide  industry  an  unusually  direct  association 
and  utility  with  programs  of  their  interest,  such  as  DAMN. 

3  3  Description  of  DF  Routines 

Concomp  personnel  also  developed  a  comprehensive 
set  of  FORTRAN  -  ca  1 1  ab  1  e  subroutines,  named  the  DF  routines^ 
(Display  File  routines),  which,  when  executed  on  the  central 
computer,  construct  DEC  338  code  to  draw  lines,  display 
labels,  and  recognize  figures  through  light  pen  detection. 
After  construction,  the  code  is  transmitted  through  a  201 
data  set  to  the  338,  which  displays  the  file  until  a  change 
is  initiated  by  the  user's  control  program.  In  constructing 
a  display  file,  control  program  tasks  include  (1)  creating 
a  construction  buffer,  (2)  constructing  figures  by  use  of 
lines  and  labels,  and  (3)  transmitting  the  file  with  an 
appropriate  name 

Other  applications  of  the  DF  routines  on  the  graphic 
terminal  include  construction  of  surfaces,  design  of  electric 
circuits,  and  simulation  of  arterial  fibrilation.  The  DF 
routines  are  also  used  in  a  number  of  computer-oriented 
courses  at  The  University  of  Michigan.  Table  3  gives  a 
brief  description  of  the  DF  routines  used  by  the  DAMN  program. 

3  4  DAMN  User  Graphics  Options 

The  primary  objective  of  the  DAMN  graphics  output  pro¬ 
gram  is  the  construction  and  display  of  a  schematic  diagram 


r 


TABLE  3 

Definitions  of  DF  Routines  Used  in  the  DAMN 
Graphics  Output  Program 


DFINI  Initializes  construction  buffer  in  central  com¬ 
puter.  All  code  generated  by  subsequent  calls 
of  DF  routines  will  be  placed  in  this  buffer 
(display  file). 

DFXYC  Draws  one  line  in  display  file. 

DFLAB  Constructs  one  label  in  display  file. 

DFOCT  Inserts  a  user- speci fied  four-digit  octal  number 
into  construction  buffer.  Used  by  graphics  out¬ 
put  program  to  change  intensities  of  lines. 

DFPJM  Allows  a  program  to  push-jump  or  call  on  a  dis¬ 
play  file  which  is  defined  as  a  subroutine.  In 
the  graphics  output  program,  three  symbols  are 
defined  in  this  way. 

DF201  The  subroutine  to  transmit  the  display  file  over 
the  AT5T  201  data  phone  to  the  DEC  338. 

DFSN  Erases  the  DEC  338  screen. 
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of  a  mechanical  system.  For  this  purpose  the  DAMN  user 
specifies  the  lines,  symbols,  vectors,  and  title  to  be 
drawn,  and  the  output  program  calls  on  the  appropriate 
DF  routines.  The  DAMN  user  has  the  option  of  defining 
whatever  lines,  vectors,  or  symbols  he  wishes  Other 
information  which  must  be  specified  in  the  data  are 
scale  factors  and  a  time  increment  which  determines 

when  the  next  frame  will  be  shown.  The  options  avail- 

\ 

able  to  the  graphics  user  and  the  information  required 
to  generate  the  display  are  summarized  in  Tables  4  and 
5  respectively 

3 . 5  Organization  of  Graphic'  Output  Program 

Figure  3  2  is  the  flowchart  for  the  graphics  output 
program  and  Figure  3.3  shows  example  output.  To  construct 
the  frames,  the  DAMN  graphics  output  program  uses  two  dis¬ 
play  files  Static  labels  such  as  the  title  and  headings 
for  velocity,  acceleration,  and  time  are  placed  in  the 
first  display  file.  As  an  example,  labels  are  shown  in 
Figure  3  3  in  the  upper  right  side  of  the  frame.  These 
labels  are  constructed,  transmitted,  and  displayed  on  the 
first  pass  of  the  output  program  and  are  branched  over  on 
subsequent  output  steps.  Also  constructed  in  the  first 
pass  are  the  three  symbols  Oi®j  and  A,  and  the  labels 

t 

for  the  vector  heads  The  vector  heads  are  constructed  by 
concatenating  the  EBCDIC  marker  identification  number  to 


TABLE  4 

Graphic  Options  and  Description  of  Output 
Available  to  DAMN  User 

Lines .  There  are  two  types  of  lines,  dark  and  light 
intensity.  By  convention,  dark  lines  or  a 
polygon  of  dark  lines  represent  rigid  parts, 
and  light  lines  represent  axes,  springs,  paths 
of  translation,  and  any  line  which  is  not  de¬ 
scriptive  of  a  rigid  part.  A  line  is  input 
by  specifying  its  two  endpoint  markers  and  its 
intensity. 

Symbols.  The  user  has  available  three  symbols:  a  1/4" 
circle,  a  1/4"  circle  crossed  in  the  middle, 
and  a  1/4"  triangle.  In  this  program  the  crossed 
circle  represents  a  rotation  which  contacts  a 
movable  part  to  ground,  a  circle  represents  a 
rotation  between  two  movable  parts,  and  a  tri¬ 
angle  represents  a  relative  translation.  The 
symbol  positions  are  defined  by  specifying  the 
marker  I.D.  number  for  the  marker  at  which  the 
symbol  is  fixed. 

Vectors.  Velocity  and  acceleration  vectors  of  any  marker 
relative  to  ground  can  be  displayed. 

Title.  A  title  of  up  to  48  characters  can  be  displayed. 
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TABLE  5 


Description  of  Data  Necessary  for  Generation 
of  DAMN  Graphic  Display 

Base.  The  user  must  specify  the  ground  marker 

I.D.  number  defined  in  the  analysis  data 
and  its  position  on  the  screen.  The 
ground  coordinates  of  all  other  markers 
are  referenced  to  this. 

Scale  Factors.  Scale  factors  are  required  to  predict 
“  the  maximum  dimensions  of  the  system, 

and  the  velocity  and  acceleration  vectors. 

Time  Increment.  The  interval  of  time  between  display  of 

~™”~*  each  new  frame  must  be  specified. 
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figure  3.2.  Mow  chart  of  the 
graphics  output  program.  Program 
i  1  o w  o ii  the  left  represents  the  con¬ 
struction  and  transmission  of  stat¬ 
ic  labels.  I’  h  e  s  e  labels  are  con¬ 
structed  at  the  first  output  step 
and  are  branched  over  on  subsequent 
output  steps.  Program  flow  on  the 
right  represents  the  construction 
o  t  each  f  r a  m  e  . 
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Figure  3.3.  Example  graphics  output.  Rigid  parts  are  repre¬ 
sented  by  straight  lines  or  polygons.  Any  number  of  velocity 
and  acceleration  vectors  can  be  shown.  In  this  example  two 
velocity  vectors  (symbols  V3  and  V5)  and  two  acceleration 
vectors  (symbols  A b  and  A7)  are  shown.  Generators  and  fields 
are  not  displayed. 


the  EBCDIC  letters  A  and  V  and  placing  the  result  in  an 
array . 

The  second  display  file  contains  the  lines  which 
describe  the  parts  in  the  system,  symbols,  and  velocity 
and  acceleration  vectors  including  their  respective 
magnitudes.  The  position  of  the  lines  (parts)  and  the 
direction  and  magnitude  of  vectors  change  with  each  frame 
and  must  be  redrawn  after  each  output  from  the  main  DAMN 
program.  The  output  program  constructs  the  system  by  re¬ 
cognizing  key  markers,  such  as  endpoints  of  lines,  reloca 
ting  these  markers  relative  to  a  base  ground  marker,  and 
scaling  all  markers  down  to  screen  size.  The  output  pro¬ 
gram  then  calls  appropriate  DF  routines  to  draw  the  speci 
fied  lines  and  symbols.  As  an  example.  Figure  3.3  is  a 
schematic  of  a  compound  pendulum  of  two  links.  Two  lines 
were  drawn  to  describe  the  links  and  two  symbols  were 
drawn  to  represent  rotations  between  one  link  and  ground 
and  between  the  two  links,  respectively. 

The  output  program  must  also  construct  velocity  and 
acceleration  vectors  if  any  are  requested.  Examples  of 
velocity  and  acceleration  vectors  are  shown  in  Figure  3.3 
Again,  the  output  program  must  recognize,  relocate,  and 
scale  the  motions  of  key  markers.  The  magnitude  and 
direction  of  the  vectors  are  computed  in  the  main  DAMN 
program,  and  the  output  program  has  only  to  scale  these 
vectors  and  call  on  the  appropriate  DF  routines. 
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Additionally,  the  output  program  must  access  the  EBCDIC 
label  f„r  vector  heads  and  call  a  DF  routine  to  add  the. 
to  the  display  Magnitudes  of  vectors  must  also  be  con¬ 
verted  fro.  floating  point  nu.bers  to  their  EBCDIC  equi¬ 
valent  and  tabulated  on  the  upper  right-hand  side  of  the 
screen  Finally  the  second  display  file  is  trans.itted 
to  complete  the  process  The  output  program  then  returns 

to  the  analysis  program  to  await  the  coordinates  for  the 
next  output  step. 

The  timing  and  quality  of  the  display  are  highly 
dependent  on  the  number  of  options  the  user  chooses,  the 
number  of  people  who  are  timesharing  simultaneously,  and 
the  number  of  dependent  variables  that  the  analysis  pro¬ 
gram  must  solve  Average  timing  experience  is  cited 
in  the  examples  which  follow 

3  6  P^cription  and  Examples,  of  Output 

Example  1  Two  degrees  of  freedom,  unconstrained,  dynamic 
compound  pendulum 

F-gure  3  4a  shows  a  sketch  of  a  compound  pendulum 
with  two  degrees  of  freedom  Each  link  has  mass  and 
moment  of  inertia  The  schematic  representation  of  this 
problem  is  shown  in  Figure  3  4b.  Parts  2  and  3  initially 
have  tero  velocity  and  are  acted  upon  by  an  unbalanced 

gravitational  force  Fivurp  •?  c  i  t 

r i gure  3  5  is  a  series  of  photo- 

graphs  of  the  compound  pendulum  as  it  is  seen  on  the 
DEC  338  Labels  V3  and  V5  define  the  magnitudes  and 


i 
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I'igure  5.4a.  Th  rec  -  d  i  mens  i  ona  1  view  of  a  two-link  compound 
pendulum  moving  in  a  gravitational  field.  Part  1  is  ground 
and  parts  2  and  5  the  two  links.  Both  parts  2  and  3  have 
appreciable  mass  and  moment  of  inertia,  and  undergo  arbitrar¬ 
ily  large  displacements.  This  is  a  two-freedom  unconstrained 
d v n a mic  svstem. 
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directions  of  the  velocities  of  the  ends  of  the  links 
Vectors  A6  and  A7  are  the  acceleration  vectors  of  the 
mass  centers  of  the  links  For  th  s  problem  DAMN  must 
integrate  a  system  of  two  second-order  differential 
equations  ’  o  obtain  the  velocity  and  position  of  each 
link  The  DAMN  output  program  draws  two  lines  and  two 
symbols,  selves  for  four  vectors  and  their  magnitudes, 
and  di-plays  3  7  frames  per  minute  a'  approximately  25% 

MTS  system  load  capacity 

Example  2  Automobile  hood  linkage 

Figure  3  6a  is  a  drawing  of  a  simplified  automobile 
hood  linkage,  and  a  schematic  of  this  linkage  i‘  shown 
m  Figure  3  6b  This  linkage  consists  of  five  moving 
parts  and  ground  (part  11  Parts  2-6  aM  have  mass  and 
mass  moment  of  inertia  A  spring  is  stretched  between 
parts  5  and  6,  with  sufficient  sTrength  to  raise  the 
hood  Also  a  torsional  spring  and  a  viscous  damper  arc 
placed  between  part-  1  and  7 

Figure  3  7  is  a  series  of  pictures  taken  at  regular 
time  intervals  during  the  motion  of  this  system  As  shown, 
the  oscillation  decays  rapidly  due  to  the  damping  Only 
one  vector,  the  velocity  vector  VI7,  is  shown,  representing 
the  velocity  ct  *  h  e  center  of  mass  situated  at  the  center 
of  the  hood  'only  the  rear  half  of  the  hood  is  shown)  A 
system  of  nine  second-order  differentia)  equations  is  solved 
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Figure  3. 6 ci.  Automobile  hood  mechanism  used  in  connecting  the 
hood  to  the  body.  The  mechanism  consists  of  six  parts.  Part 
1  is  the  automobile  body  and  part  6  the  hood.  A  spring  is 
stretched  across  parts  5  and  6  and  is  sufficiently  strong  to 
raise  the  hood  from  a  near-horizontal  position.  The  motion  of 
the  hood  is  impeded  by  a  rotational  viscous  damper  (not  shown) 
between  part  1  and  part  2.  This  is  a  single-freedom,  con¬ 
strained  dynamic  system. 


Figure  3.6b.  Two-dimensional  schematic  representation  of  the 
hood  mechanism  of  Figure  3.6a.  Rotations  (pin  joints)  are 
represented  by  small  circles;  parts  are  represented  by  straigh 
lines  and  polygons. 


nc remen ted  in 
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as  described  in  Section  2  and  the  output  program  dis¬ 
plays  7  lines,  7  symbols,  and  the  velocity  vector. 

Other  Examples 

Figures  3.8a  through  3.8d  show  other  systems 
which  have  been  simulated  by  DAMN.  These  examples 
show  that  DAMN  can  accept  a  wide  variety  of  situations 
such  as  rotational  and  translational  contact,  kinematic 
constraint,  multifreedom  dynamic  systems,  constrained 
dynamic  problems,  inertially  driven  problems,  etc. 

These  and  other  examples  also  showed  that  the  number 
of  frames/min.  was  more  often  dependent  on  data  trans¬ 
mission  time  than  an  problem  integration  time.  Exper¬ 
ience  has  also  shown  that  graphics  is  particularly  use¬ 
ful  in  analy2ing  systems  which  undergo  large  deflections 
and  is  less  interesting  for  small-scale  vibration. 

3  7  Man-Machine  Interaction 

An  additional  feature  of  the  DAMN  graphics  output 
program  is  man-machine  interaction  through  the  use  of 
an  interpretive  command  language  This  is  not  yet  fully 
developed,  but  at  present  a  user  can  interrupt  the  dis¬ 
play  at  any  point  and,  by  issuing  an  appropriate  command, 
terminate  the  program,  access  a  new  data  set,  alter  the 
display,  or  change  system  parameters.  When  the  user 
wishes  to  modify  the  display  he  first  issues  additional 
subcommands  Through  these  subcommands  he  adds  or  deletes 


Figure  3.8.  Photographs  of  the  DEC  338  graphic  display  screen 
showing  other  systems  simulated  by  DAMN. 


Figure  5.8a.  A  kinematic  two-slider  mechanism  driven  by  a 
constant  velocity  input  motion.  Vectors  V 7  and  V12  are  the 
velocities  of  the  sliders. 


Figure  3.8b.  \  one- freedom  four -bar  linkage  (no  input  mo¬ 

tion).  Each  link  lias  mass,  moment  of  inertia,  and  initial 
velocity.  Vectors  \  1 0  and  A 10  are  the  velocity  and  accelera¬ 
tion  of  the  mass  center  of  the  coupler. 


Figure  3.8c.  A  simplified  two-dimensional  automobile  suspen¬ 
sion  system.  The  canted  line  is  the  frame,  supported  verti¬ 
cally  by  two  springs  (not  shown),  one  at  each  end.  The  ver¬ 
tical  line  shows  vertical  displacement  of  the  center  of  mass, 
and  provides  a  reference  for  a  rotational  damper  (not  shown). 


Figure  3.8d.  A  two-dimensional  simulation  of  a  washing 
machine.  The  rotating  tub  is  represented  as  a  sinusoidal  os¬ 
cillator  mounted  on  a  base-plate  suspended  from  the  frame  by 
links.  The  effect  of  an  eccentric  mass  in  the  tub  is  repre¬ 
sented  by  specifying  a  mass  for  the  oscillator. 
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lines,  symbols,  or  vectors,  alters  scale  parameters, 
and  changes  the  frame  increment.  This  feature  allows 
the  user  to  add  a  temporarily  significant  vector  for 
a  particular  interval  of  time  and  then  delete  it  to 
simplify  the  display.  As  planned,  the  command  language 
will  allow  the  user  to  add  or  delete  fields  or  gener¬ 
ators,  and  to  alter  part  geometry  and  field  or  generator 
parameters,  then  restart  the  system,  and  examine  the 
effects  of  the  change. 
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